Molecular spinning by a chiral train of short laser pulses 
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We provide a detailed theoretical analysis of molecular rotational excitation by a chiral pulse train 
- a sequence of linearly polarised pulses with the polarisation direction rotating from pulse to pulse 
by a controllable angle. Molecular rotation with a preferential rotational sense (clockwise or counter- 
clockwise) can be excited by this scheme. We show that the directionality of the rotation is caused by 
quantum interference of different excitation pathways. The chiral pulse train is capable of selective 
excitation of molecular isotopologues and nuclear spin isomers in a mixture. We demonstrate this 
using 14 N2 and 15 N2 as examples for isotopologues, and para- and ortho-nitrogen as examples for 
nuclear spin isomers. 

PACS numbers: 32.80.Qk, 33.80-b, 37.10.Vz, 42.65.Re 



I. INTRODUCTION 

The control of rotational molecular dynamics by non- 
resonant strong laser fields has proven to be a powerful 
tool. It allows creating ensembles of aligned [lmll , ori- 
ented [H, 0| or planarly confined molecules [sl-tl3jT The 
proposed and realised applications are numerous, includ- 
ing control of chemical reactions @, [lij] , high harmonic 
generation [ID, [l|| , control of molecular collisions with 
atoms [13 or surfaces [Hill, and deflection [H-ii] of 
molecules by external fields. 

An important challenge for strong field rotational con- 
trol is selective excitation in a mixture of different molec- 
ular species. Isotopologue selective control was demon- 
strated in [26M28| , using constructive and destructive in- 
terference induced by a pair of delayed laser pulses [111 . 
With a similar scheme also nuclear spin isomer selec- 
tive excitation was achieved [U E( • More recently, 
isotopologue selective rotational excitation by periodic 
pulse trains has been demonstrated [HJ and connection 
of this scheme to the problem of Anderson localisation 
was shown [33| . 

One current direction of strong field rotational control 
focuses on the excitation of molecular rotation with a pre- 
ferred sense of the rotation. In the "optical centrifuge" 
approach [H, |H|-[36|, the molecules are subject to two 
counter-rotating circularly polarised fields, which are lin- 
early chirped with respect to each other. The resulting 
interaction potential creates an accelerated rotating trap, 
bringing the molecules to a fast spinning state. The al- 
ternative "double pulse" scheme reaches the same go al by 
using two properly timed linearly polarised pulses [9|, llCf . 
The first pulse induces molecular alignment. When the 
alignment reaches its peak, the second pulse, whose po- 
larisation is rotated by 45 degree with respect to the first 
one, is applied and induces the directed rotation. More 
recently, an approach using a "chiral pulse train" was 
demonstrated [371 ] . in which a train of pulses is used, 
where the pulse polarisation is changed by a constant 
angle from pulse to pulse (see Fig. [1}. 

In this article, we provide a detailed theoretical anal- 
ysis of the rotational excitation by a "chiral pulse train" 



demonstrated in |37j . A thorough description of the ex- 
perimental procedure is presented in a companion arti- 
cle [38| . 

For the present article, the structure is as follows. In 
Section ID we introduce the model for the lascr-molccule- 
interaction. Then, we consider excitation scenarios for 
two kinds of molecules. The first one is N2, representing 
a simple diatomic molecule which is well described by the 
standard model of a rigid rotor. The second molecule we 
consider is O2. Unlike the nitrogen molecule, it has a 
non-zero electronic spin in its ground state, leading to 
a more complex structure of the rotational levels. Our 
analytical and numerical results are presented in Sec- 
tion UTTJ Here, we first show the results for the excitation 
of 14 N 2 by a chiral train of equally strong pulses. Next, 
we demonstrate the prospects of selective excitation of 
nuclear spin isomers and isotopologues by such trains. 
Finally, the results for oxygen molecules interacting with 
a train of non-equal pulses are shown and compared with 
the experiment [37| . In the last section, we summarise 
the results and conclude. 



II. MODEL AND NUMERICAL TREATMENT 
A. Model 

We consider the following scenario: A train of ultra- 
short laser pulses interacts with a gas sample of linear 
molecules like N2 or O2. The pulses are applied with 
a constant delay t between them. Each pulse is lin- 
early polarised, but the polarisation vector is rotated 
from one pulse to the next one by the angle S, such that 
the whole pulse train rotates with a rotational period 
Ttrain = 2ttt / S (see Fig. []}. Choosing the laser propa- 
gation axis as z-axis, the electric field of the n th pulse is 
given as 



E n (t) = £ n (t)e n cos(w£ 



(1) 



where e„ = (cosn<5, sinnJ, 0) is the polarisation vector, 
uj is the carrier frequency, and 4> is the phase. We con- 
sider a laser for which the carrier frequency is far detuned 
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FIG. 1. (Colour online) Sketch of the considered scenario: 
A train of linearly polarised laser pulses interacts with linear 
molecules. The polarisation axis is rotated by an angle S 
between each pulse, and the time-delay r between the pulses 
is constant. 



from any electronic or vibrational resonance. The laser 
pulses therefore interact with the molecules via Raman- 
type excitations of the rotational levels. The envelope of 
the electric field is given as 



£ n (t) = E n cxp [-(* - tit) 2 /(2a 2 )] 



(2) 



Here, a determines the pulse duration. 

The non-resonant laser pulse induces a dipole in the 
molecule via its electric polarisability, and then interacts 
with this induced dipole. Averaging over the fast os- 
cillations of the electric field, we arrive at the effective 
interaction potential [39| 



V 



4 (Aacos 2 /?,, 



(3) 



Here, Aa = a« — a± is the polarisability anisotropy of 
the molecule, where aii and a± are the polarisabilities 
along and perpendicular to the molecular axis, respec- 
tively. The angle j3 n is the angle between the molecular 
axis and the polarisation direction of the n th pulse. The 
last term in Eq. ([3]), is independent from the molec- 
ular orientation and does not influence the rotational dy- 
namics. We will therefore omit it in the following. 

It is convenient to introduce an effective pulse strength 
P, which corresponds to the typical amount of angular 
momentum (in units of K) transferred to the molecule by 
the pulse. For a single pulse, it is defined as 



(4) 



Here, I pea k is the peak intensity of the pulse, c is the 
speed of light, and e is the vacuum permittivity. 

In this work, we consider two kinds of pulse trains. 
The first one is a train of N equally strong pulses, such 
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FIG. 2. Effective interaction strengths for a train with modu- 
lated intensities as described by Eq. ((6}. The total interaction 
strength is P to t = 1. 



that the effective interaction strength P n of the n th pulse 
is given as 



Pn 



Ptot/N, 



(5) 



where P tot = P n is the total strength of the whole 
pulse train. Such a pulse train can be generated e.g. 
by nested interferometers EH • The second kind are 
trains like the ones used in the experiments [13, [H[, 
which were created by pulse shaping techniques. In this 
case, the effective interaction strength of the n th pulse is 
given as 



Pn — PtotJn(A) 



(G) 



where J n is the Bessel function of the first kind, and A 
is a parameter. Since J n (x) -C 1 for \n\ > \x\, this train 
contains about 2A+1 non-zero pulses. In Fig.[2]we depict 
the intensity envelope of the train for different values of 
A. 



B. Numerical treatment 

1. Nitrogen 

At first we consider molecular nitrogen as an example 
of a simple linear molecule. Since the laser pulses are 
assumed to be far off-resonant from electronic or vibra- 
tional transitions, it is sufficient to consider only the rota- 
tional excitation in the vibronic ground state. The rota- 
tional eigenfunctions are the spherical harmonics |J, M). 
Here, J is the total angular momentum, and M is its pro- 
jection on the Z-axis, which we have chosen to be along 
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the laser propagation direction. Note that for N2 in its 
electronic ground state, the total angular momentum J is 
equal to the orbital angular momentum O of the rotation 
of the nuclei. Although we are interested in the latter, for 
simplicity we keep to the more common notation using 
the total angular momentum J. The rotational levels are 
given as Ej = BJ(J + 1) - DJ 2 (J + l) 2 , where B is the 
rotational constant and D is the centrifugal distortion 
constant. 

For the numerical treatment of the problem it is conve- 
nient to express the wave function as a linear combination 
of the rotational eigenfunctions: 

!*(*)> = ^ C ^)e- iEjt/h \ J, M). (7) 

J,M 

Inserting the expansion ([7]) and the interaction poten- 
tial ([3]) into the time-dependent Schrodinger equation, 



Here, we use the relations 



«« - mm 



(8) 



we obtain 

ih J2 9Cj '^' {t) e- iE ^ h \J',M') 



J'M' 



C JIMI (t)e- iE ^/ h V(t)\J',M'). (9) 



J'M' 



Multiplying from the left by ^-(J, M\e iEjt/h we obtain 
a set of coupled differential equations for the expansion 
coefficients Cj t M(t): 

dC JM (t) 
dt 

=lr Y C, rM '{t)e-^ E -''- E ^ t ' h {J,M\V{t)\J',M') 



J',M' 
+00 



ih 
Act 

n= — oo J',M' 

x (J, M\ cos 2 p n \J',M') 



J2 £2(t) E Cj,Mt)e-« B "- BjWn 



(10) 



Here, j3 n is the angle between the molecular axis and the 
polarisation direction of the n th pulse. 

The matrix element (J,M\ cos 2 (3 n \ J', M') is obtained 
as follows. First, cos 2 (3 n is expressed as 

cos 2 (3 n = cos 2 (n<5) sin 2 9 cos 2 4> + sin 2 (nS) sin 2 9 sin 2 <j> 
+ i sin(2n<f) sin 2 9 sin(2<?!>) , (11) 

where 9 and <f> are the polar and azimuthal angle of the 
molecular axis, respectively. Then, we express cos 2 /3 n in 
terms of the Wigner rotation matrices [42[ O^lv as 



cos 2 /? - ---D {2) * + —e t2nS D {2) * + —e- 
cos Pn — 2 3 00 y/Q \/6 



i2nS r>( 2 )* 
-^20 • 

(12) 



sin 2 9 cos 2 1 



sin 2 9 sin 2 1 



3s 



D (2) * 
V6 V 20 



v / 6 V 20 



(2)* 
-20 



2)* 
-21) 



-D, 



(2)* 

no 



3 00 + 3 



sin 2 9 sin(20) = 
Finally, by using (42j 



(2)* 



- D 



(2). 



(13a) 



(13b) 



(13c) 



{J,M\D^*\J',M') = (-1) M V(2JT 
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J o')( 



1)(2J' 

J 2 J' 

■M Af M' 



(14) 



where the brackets denote the Wigner 3-j symbol, we 
obtain the matrix element ( J, M\ cos 2 fi\ J', M'). Note 
that only levels with A J = 0, ±2, ±4, ... and AM = 
0, ±2, ±4, ... are coupled. 

In our simulations, we solve Eq. (|10p numerically. We 
do ensemble averaging by solving Eq. (fT0|) for different 
initial states {^initial) = I^Oj-^o) and weighting the re- 
sult by the Boltzmann factor of the initial state. Note 
that the Boltzmann factor includes a degeneracy factor 
arising from nuclear spin statistics [43|. For example, 
the nitrogen isotope 15 N has a nuclear spin of I = 1/2. 
Therefore, the diatomic molecule 15 N2 can have a total 
nuclear spin of J = 1 (ortho-nitrogen) or / = (para- 
nitrogen). The former has three degenerate nuclear spin 
wave functions, which are symmetric with respect to an 
exchange of the two nuclei, and the later has one antisym- 
metric nuclear spin wave function. Due to the fermionic 
nature of 15 N, the total wave function of the molecule 
has to be antisymmetric with respect to the exchange of 
the nuclei. Therefore, ortho- and para-nitrogen can be 
distinguished by their rotational wave functions: Ortho- 
nitrogen is only found with odd angular momentum J, 
para-nitrogen only with even angular momentum J, and 
the ratio of even to odd states is 1:3 due to the degener- 
acy of the nuclear spin wave functions of ortho-nitrogen. 
For 14 N with a nuclear spin of / = 1, there are three nu- 
clear spin isomers, two with symmetric nuclear spin wave 
functions (one of them five- fold degenerate) , and one with 
three-fold degenerate antisymmetric nuclear spin wave 
functions. The resulting ratio of even to odd rotational 
states is 2:1. 



2. Oxygen 

Molecular oxygen has a more complex rotational spec- 
trum than simple diatomic molecules like nitrogen or hy- 
drogen [13, Sj|. The electronic ground state is a 3 E~ 
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FIG. 3. The lowest rotational levels of O2 in its electronic 
and vibrational ground state. The energies are given in cm -1 . 
Also shown are the allowed transitions between these levels 
induced by the laser pulse (J3J). 



state, so the total electronic spin S is non-zero. This 
gives rise to spin-spin and spin-orbit coupling, and there- 
fore the total angular momentum J docs not solely de- 
scribe the nuclear rotational motion as for N2, which has 
a 1 S+ electronic ground state. In particular, J = N + S, 
where N is the orbital angular momentum; since the elec- 
tronic orbital angular momentum is zero, N is identical 
to O, the nuclear orbital angular momentum. The ro- 
tational quantum number J can take the values J = 
N + S,N + S - 1,...,N - S. Therefore, for oxygen in 
its vibronic ground state with S = 1, each level N is 
split into three levels with J = N — 1, N, N + 1, as is 
shown in Fig. [3] The splitting is stronger for low val- 
ues of TV. Additionally, for symmetry reasons, only odd 
values are allowed for N 

For the numerical treatment, we express the wave func- 
tion as a linear combination of Hund's case (b) basis 
states El El: 



E 

JNM 



CjNM^e-^-'^/^K; NA; NSJM) . (15) 



Here, A is the projection of the electronic angular mo- 
mentum on the molecular axis, N is the orbital angular 
momentum, S is the electronic spin, J is the total angu- 
lar momentum, M is the projection of the total angular 
momentum on the Z-axis, and rj is a combined quan- 
tum number of the remaining vibronic quantum numbers. 
Ejn are the energies of the rotational states, see Fig. [3] 
As before, we assume that the molecules are initially in 
their vibronic ground state. Since the interaction does 
not induce any vibronic transitions, it is independent of 
r), and furthermore A = and S = 1 are constant. For 



ease of reading, in the following we denote the eigenstatcs 
in short as \t]A; N A; NSJM) = \ JNM). 

As before, we insert the expanded wave function (|15p 
into the time-dependent Schrodinger Equation and ob- 
tain a system of differential equations for the expansion 
cocffccients Cjnm- 



dC. 



JNM 



(*) 



at 



Aa \ -> „ , . 
J 7F E £ n® 

C.,, N , M >{t)e- t{E ->'»- EjN)t/h 
x (J, N, M\ cos 2 P n \J', N', M') . 



+ OO 



n— — 00 

E 

J'N'M 1 



(16) 

The matrix elements ( JNM\ cos 2 f}\ J'N'M') are deter- 
mined as follows. Firstly, we use Eq. (fT2| to replace 
cos 2 j3, which yields 

(J N M\ cos 2 fi\ J'N'M') = 

-{JNM\J'N'M') 



--{JNM\D$*\J'N'M' 



(17) 



+J-e i2nS (JNM\D l ?l*\ J'N'M') 
V6 

+ -L e - l2nS (JNM\D<$*\ J'N'M') . 
V6 

We now have to determine the value of 
(ANSJM\D M l* \AN'SJ'Al'). Here we explicitly 
write down all quantum numbers (apart from rf). We 
use the Wigner-Eckart-theorem (see e.g. [42[) to exclude 
the dependence on the molecular orientation: 



(ANSJM\D M } * Q \AN'SJ'M') = (- 



,J-M 



J 



J' 



-M M M' 



(ANSJ\\Df*\\AN'SJ'), (18) 



where the round brackets is the Wigner 3-j symbol. The 

dot in the subscript of the rotation matrix indicates that 

this matrix element is reduced regarding the orientation 

in the space-fixed coordinate system j42(. Next, we use 
(2)* 

the fact that Dq' does not act on the electronic spin, 
so we can exclude S from the matrix element as well and 
obtain EH 

(ANSJ\\D { q ] *\\AN'SJ') = 

(-l) J ' +N+2+s ^(2J + 1)(2 J' + 1) 

N j J N S 2 yAN\\Df*\\AN') . (19) 

Here, the curly brackets denote the Wigner 6-j symbol. 
Finally, the reduced matrix element in Eq. (|19p is given 
as 



(AiV||^ ( 2) *||AiV') = 
(~1) N - A y/(2N + 1)(2N' 



N 2 N' 
-A A 



(20) 
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Here, we used Eq. (5.186) in [42| and applied it for Hund's 
case (b). We insert now 5 = 1 and A = 0, and obtain 
the matrix elements of the rotation matrices as 

(ANSJM\D$* \AN'SJ'AI') = 

x y/(2J + 1)(2 J' + 1){2N + 1)(2N' + 1) 



(21) 



J 


2 


J ' ) 




2 


N' \ 


-M 


M 


M J 


o 





j 



N' J' 1 
J N 2 



(-1) 



J+J'-M+l 



Inserting Eq. (f2~Tj) into (fTT|) yields the matrix elements 
(JNM\ cos 2 /3| J'N'M 1 ). In order to lower the numerical 
effort, we treat the pulses as delta-pulses (sudden approx- 
imation), i.e. we neglect the molecular rotation during 
each pulse. Comparison with experiment (37j shows that 
this approximation is well justified for pulses of a dura- 
tion of 500 fs. Using the method of an artificial time 
parameter £ as described in Q , the differential equations 
for the expansion coefficients for a single laser pulse be- 
come 



dC JNM {Q 



J',N'M' 

x (J, N,M\ cos 2 p n \j',N',M') , (22) 



where P n is the effective interaction strength introduced 
above. Setting Cjnm{£, = 0) to the values just before the 
pulse, we obtain the expansion coefficients right after the 
pulse as Cjnm{£, = 1) |2|- To obtain the final expansion 
coefficients after the whole pulse train, we solve Eq. ([22]) 
for every pulse, letting the wave packet (|15j) evolve freely 
between the pulses. To account for thermal effects, we 
do ensemble averaging over the initial state. Since 16 
has a nuclear spin of / = 0, there are no degeneracies due 
to the nuclear spin wave functions. However, only odd 
values are allowed for the orbital angular momentum N . 



III. RESULTS 

We will first present the results for excitation of nitro- 
gen molecules by a train of equally strong pulses. We will 
then demonstrate how such pulse trains can be used to 
selectively excite isotopologues and nuclear spin isomers 
in molecular mixtures. Finally, we will show results for 
the excitation of the more complex oxygen molecules by 
a train of unequal pulses (given by Eq. ([6|)), in order to 
compare our results with recent experiments (37l . l38j . 

We define the final population Q(J) of a rotational 
level J as 



JM 



(23) 



A I 



Here, i denotes the initial state and gi is its statistical 
weight. We also define the directionality e(J) of the ex- 



cited wave packet as 



e(J) 



Ql(J)-Qb(J) 
Ql(J) + Qr(J) 



(24) 



where Ql(J) and Qr(J) are the counter-clockwise rotat- 
ing and the clockwise rotating fraction of the population 
of the level J, 



Ql(J) = X>( £ \Q 

i \M>0 
i \M<0 



m| 2 + 1/2|C v/0 | 2 J (25a) 
r ^M| 2 + l/2|a,Jo| 2 ) ■ (25b) 



Note that the states with M = are accounted half for 
clockwise and half for counter-clockwise rotation. A pos- 
itive (negative) e(J) indicates a preferentially counter- 
clockwise (clockwise) rotation. 



A. Excitation of nitrogen molecules with a train of 
equally strong pulses 

In the following, we show the results for 14 N2 molecules 
interacting with a train of eight equally strong pulses 
with durations of a = 30 fs (see Eq. ([2])) and a total 
interaction strength of Ptot — 5. The peak intensity of 
a single pulse is therefore approximately 3 • 10 12 W/cm 2 . 
The pulse duration is well below the rotational periods of 
the highest expected excitations (remember that P cor- 
responds to the typical angular momentum in the units 
of H transferred by the pulse). The molecules are con- 
sidered being initially at a temperature of T = 8 K. At 
this temperature there is a considerable initial (thermal) 
population in the level J = 2, with Qth{2) = 0.25. Also 
in J = 3 there is some initial population, Qt^(3) = 0.02. 
The levels J = 4 and J = 5 are not populated (note that 
due to nuclear spin statistics, two thirds of the popula- 
tion is found in the even levels, and one third in the odd 
ones). 

In Figure [4] the population Q(J) is shown for the ro- 
tational levels J = 2, 3,4, and 5 for 14 N2 molecules, and 
Figure [5] displays the directionality e(J) for the same lev- 
els. The plots show the population and the directionality 
as a function of the pulse train period r and the pulsc- 
to-pulse polarisation angle shift S. 

The population plots for all levels show a distinct pat- 
tern of diagonal and horizontal lines. These lines are 
described by the equation 



texc(J) m + AM 



2- 



(26) 



Here, m is an integer and AM = 0, ±2 (0 yields the hor- 
izontal lines, +2 corresponds to the diagonal lines with 
a positive slope, and —2 yields the diagonal lines with 
a negative slope). Furthermore, t exc is the period cor- 
responding to the excitation from the level J — 2 to the 
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FIG. 4. (Colour online) Population Q(J) (see Eq. ((23j)) of 
different rotational levels J for 14 N2 at T = 8K, after inter- 
acting with a train of eight equal pulses. The total interaction 
strength is Ptot = 5 and the pulse duration is a — 30 fs. The 
minimum of the colour-coding is for all panels, and the 
maximum is 0.6 for (a), 0.25 for (b) and (c), and 0.1 for (d). 



FIG. 5. (Colour online) Directionality e(J) (see Eq. (J24j|) of 
different rotational levels J for 14 N2 at T = 8K, after inter- 
acting with a train of eight equal pulses. The total interaction 
strength is Ptot = 5 and the pulse duration is a = 30 fs. Pos- 
itive e corresponds to counter-clockwise rotation, negative e 
corresponds to clockwise rotation. 



level J, and is given as 

t exc (J) = 2nh/{Ej - Ej_ 2 ) = t rev /(2J - 1) , (27) 

where t rev = tin/ B is the rotational revival time (8.38 ps 
for 14 N2 in its vibronic ground state). 

The directionality plots in Fig. [5] show in general the 
same structure as the population plots, although now 
the horizontal lines are missing. Furthermore, we can see 
that the diagonals with a positive slope correspond to 
a counter-clockwise rotational sense (e(J) > 0), and the 
diagonals with a negative slope to a clockwise rotation. 
Also, next to the main diagonals, there is a chessboard 



pattern visible, especially for higher levels. 

The general structure of the population and direction- 
ality plots is ver y s imilar to experimental observations 
(Figs. 5 and 6 in [38|), although the latter do not resolve 
the fine chessboard pattern. It should be noted that in 
the experiments [38[ a train of unequal pulses described 
by (O was used, whereas the results presented here are 
for a train of equally strong pulses. 

The structures seen in Figures [4] and [5] can be explained 
as the result of the quantum interference of different exci- 
tation pathways, as we will show now. For simplicity, we 
will treat the pulses as delta-pulses in the following anal- 
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ysis. This is well justified, as the utilised pulse duration 
of a = 30 fs is much shorter than the relevant rotational 
periods of 14 N2 . The evolution of the wave packet over 
one period of the pulse train is given by 



state to the final state after the last pulse as 

N 

U = JJ e - mSj * ,h We inS3jh e- 1 & T . (30) 



|*(f+)) = e ^co S 2 /3„ e - 4 Pr/(2/ft)|^^+_ i ^ 



(28) 



where i+ is the time-instant right after the n th pulse, J 
is the angular momentum operator, and / is the moment 
of inertia. The interaction term can be expressed as 



iP cos [3 n 



R(n5, 7r/2,0)e 



iPc 



'ie _x (0,-7r/2,-n5) 



iSJg/H -ij v w/(2h) iP cos 2 ij y ir/(2h) inSJ^/h 



(29) 



Here, Br 1 rotates the basis from the space fixed sys- 
tem (quantisation axis along the laser propagation) to a 
"pulse fixed" system (quantisation axis along the electric 
field polarisation of the n th pulse). The operator W is 
the same for every pulse. Using (f2"§)) . we can express the 
evolution operator that brings the system from its initial 



The probability of the transition from | J'M') to \ JM) is 

2 



given as 



(JM\U\J'M') 



Using the expansion 



JP cos 2 = 1 + i p cos 2 0_;P og 4 ( 



we can express the evolution operator as 

N 

u=n e -^A e -*-w(2fc) (x + zPcoS 2 6 + ) 



(31) 



n=l 



x e 



i.J y Tr/(2h) in5J z /h -i^T 



rN 



N 



-i-^jjrr{N — n) —in8J z /h 



+ iP^^e -2i»-\" '-'e 

-ij y n/(2h) CQg 2 g e iJ y n/(2h) ^inb J z I fi e -i4m Tn 

(32) 



x e 



The approximation in the last line is valid in the limit of 
weak pulses (P < 1). In the following, we will only con- 
sider this limit. Using Eq. (|3"2"|) as evolution operator, the 
total probability for a transition from the state | J'M'} to 
another state | JM) is given as 



(JM\U\J'M r 



p2 



N 



e -i(N-n)EjT/h^j M |y| j' M ^ e -i(M-M')n6 e -inE,,T/h 



n=l 



N 



= P 2 \(JM 



V 



J'M')\ 2 J2 cos 



n.n' — l 



Y AEt 



AMS) (n - n') 



(33) 



Here, V = e-^v^/C 2 ^ cos 2 ^^^ 2 ^ and AE = Ej - 
Eji . The term <£■, and therefore the transition amplitude, 
is maximised, if the first factor in the argument of the 
cosine is a multiple of 27r, which yields 

where m is an integer. This condition is equivalent to (|26p 
and exactly describes the lines in Figures S] and O 

Using these insights, we can now explain the results 
seen in Fig. 2] and [51 The patterns are the result of quan- 
tum interferences of different excitation pathways. These 
interferences are constructive, when the condition (|34|) 



is fulfilled, causing the lines seen in Fig. [4] and By 
the help of Eq. (j34|) we can also see that the horizontal 
lines are due to transitions with no change of the pro- 
jection M, i.e. AM = 0. There are no horizontal lines 
in the directionality plots, since AM = means that 
there is no change in the sense of the rotation. The di- 
agonals with a positive slope are due to transitions with 
AM = 2. The increase of M shifts the rotational sense 
towards a counter-clockwise direction, and therefore in- 
creases the directionality e(J). The opposite is found for 
the diagonals with a negative slope, which correspond to 
AM = -2. 

The chessboard pattern seen in the directionality plot 
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detuning x (in units of t ) 

FIG. 6. (Colour online) Sum |<J>| of the interference terms 
in Eq. (133j> as a function of the detuning x from the condi- 
tion (|34p . The plot shows <E> for three different numbers N of 
pulses, and the result is normalised to TV 2 . 



can be explained by looking at the sum $, when t is 
detuned from the condition (|34l) : 

t = t exc AM— + x . (35) 

Here, t exc = 2nh/ AE. In Figurc[6]we plot $ as a function 
of the detuning x. It can be seen that next to the main 
peaks at integer x/t exc , there are weak oscillatory beats 
in between. The minima of those beats are found at 
x/texc = m/N, where N is the number of pulses and m 
and TV are mutually prime. These "side bands" are weak 
and therefore they can not be seen in the population plots 
in Fig.[4j On the other hand, the directionality measures 
the relative difference of the populations, so these weak 
side bands become visible, if the thermal population of 
the rotational level is sufficiently small. 

From Figure [6] we can see that an increase of the num- 
ber of pulses leads not only to an increase in the number 
of "side bands" , but also to a narrowing of the main peak 
of $(x). Therefore, we expect a narrowing of the lines 
seen in the population plots for larger N. This can be 
seen in Fig. [7] for the population of the levels J = 2 and 
J = 3. Here, we use the same parameter values as in 
Fig. |U but twice as many pulses, while keeping the total 
interaction strength Ptot constant. A similar effect of the 
narrowing of the resonance when increasing the number 
of pulses was already found for the quantum resonance 
at the full rotational revival (33|. 

B. Selective Excitation 

1. Nuclear spin isomer selective excitation 

The nitrogen isotopologue 15 N 2 can be found as ortho- 
nitrogen with a total nuclear spin of I = 1, or as para- 
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pulse to pulse polarisation shift 8 (% rad) 

FIG. 7. (Colour online) Population Q(J) (see Eq. p3]l) of 
the rotational levels J = 2 and J = 3 for 14 N 2 at T = 8K, 
after interacting with a train of 16 equal pulses. The total 
interaction strength is P to t = 5 and the pulse duration is 
a = 30 fs. 



nitrogen with a total nuclear spin of / — 0. These spin 
isomers can be distinguished by their rotational wave 
functions [43] ]: Ortho- nitrogen is only found with odd 
angular momentum J, para-nitrogen only with even an- 
gular momentum J. 

In Figure [5] we show the final population for the lowest 
rotational levels, after excitation by a pulse train with 
5 = and period r close to one quarter of the revival 
time (marked by the dashed line). By tuning the time- 
delay between the pulses one can choose which state is 
excited strongest. Moreover, for t < t rev /4 only odd lev- 
els are significantly excited, whereas for r > t rev /4 only 
even levels are significantly excited. Therefore, by choos- 
ing r slightly smaller (slightly larger) t rev /4, one selec- 
tively excites ortho- nitrogen (para- nitrogen) . This effect 
is shown in Figure [9] (a) , which displays the absorbed en- 
ergy for both spin isomers. One can see separate peaks 
for both isomers, in particular close to r = t rev /A and 
t = 3t rev /A. This selective excitation of spin isomers 
around a quarter of the revival time was demonstrated 
in a recent experiment [32j . 

At 5 0, one may use chiral pulse trains to bring dif- 
ferent spin isomers to a rotation of opposite sense. The 
best selectivity is achieved for 6 = tt/4. In particular, at 
S = 7r/4 and r = t rev /4 all excited even states have a pos- 
itive directionality, and all excited odd states have a neg- 
ative directionality. The reverse is found at r = 3t rev /4. 
The opposite directionality of even and odd rotational 
states at r = t rev /4 and r = 3t re „/4 was also demon- 
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FIG. 8. (Colour online) Final population of the rotational 
levels of 15 N2 after interaction with 8 equal pulses with a 
duration of a — 30 fs and a non-rotating polarisation (8 — 0) 
and Ptot — 5, at T = 8 K. Each curve is normalised to its 
maximum. The dashed line indicates the quarter revival time 
ireu/4. For t slightly smaller than t re v/4, only odd states are 
excited, for r slightly larger, only even states. 




4 6 
time-delay x (ps) 




4 6 
time-delay x (ps) 



strated experimentally (see Fig. 7 in [HI). This effect 
allows for spin isomer selective excitation, as shown in 
Fig. [5] (b). Here, the projection of the angular momen- 
tum on the z-axis is shown for both spin isomers. For 
t = t rev /A, ortho-nitrogen exhibits counter-clockwise ro- 
tation, whereas para-nitrogen rotates clockwise. The op- 
posite is found at r = 3t rev /A. 

2. Isotopologue selective excitation 

Isotopologucs arc chemically identical molecules with a 
different isotopic composition, e.g. 14 N2 and 15 N2. Due 
to the different moments of inertia, isotopologues have 
different rotational time-scales. For example, 14 N2 has 
a rotational revival time of t rev = 8.38 ps, whereas for 
15 N2 the revival time is t rev — 8.98 ps. Using that the 
rotational excitation is strongest if the pulse train pe- 
riod r equals the rotational revival time, we can selec- 
tively excite 14 N 2 or 15 N 2 by tuning the train period to 
t = 8.38 ps or t = 8.98 ps, respectively [UHl]. In 
Figure [TU] (a) we show the absorbed rotational energy of 
the two nitrogen isotopologues after interaction with a 
non-rotating (6 = 0) pulse train, and one can clearly see 
the selective excitation at the respective revival times. 

Inducing counter-rotation of different isotopologues is 
more challenging. For spin isomers the time-scales were 
identical, and one could use the different directionality 
of even and odd states for r = t rev /A to induce counter- 
rotation. For isotopologues the time-scales are different. 
Counter-rotation can only be excited if for some set of 
parameters the pulse train accidentally excites rotation 



FIG. 9. (Colour online) Selective excitation of the nuclear 
spin isomers of 15 N2 by interaction with a periodic pulse train 
of eight equally strong pulses with a duration of a = 30 fs. 
The molecules are initially at T = 8 K. The total interaction 
strength is Ptot = 5. The results are shown as a function 
of the pulse train period r. (a) Absorbed rotational energy 
for a non-rotating (S = 0) pulse train, (b) Projection of the 
molecular angular momentum on Z for a rotating (S — 7r/4) 
pulse train. The angular momentum is given in units of h. 
The dashed vertical lines indicate quarter, half, three quarter 
and full rotational revival times. 



of opposite direction in the isotopologues. For 14 N2 and 
15 N2 three such regions can be seen in Fig. \W\ (b) (see 
shaded region): At r ss 2.2 ps and r « 7.65 ps the heav- 
ier isotopologue rotates predominantly counter-clockwise 
((J z ) > 0), and the lighter isotopologue rotates mainly 
clockwise; at tw 8.65 ps the opposite is found. 



C. Oxygen molecules in a chiral pulse train 

As a special example, we now consider the excita- 
tion of oxygen molecules by a chiral pulse train. In- 
stead of identical pulses we consider the more complex 
pulse sequence ([6]), corresponding to the one used in 
experiment [13] (see Fig. [2]). We also use parameters 
corresponding to this experiment: The total interaction 
strength is P tot = 7.5, and A = 2 (see Eq. ©). There- 
fore, the strongest pulse in the train has an effective inter- 
action strength of P = 2.5, which corresponds to a peak 
intensity of approximately 8 • 10 12 W/cm 2 . Due to the 
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FIG. 10. (Colour online) Selective excitation of the nitrogen 
isotopologues 14 N2 and 15 N2 by interaction with a periodic 
pulse train of eight equally strong pulses with a duration of 
a = 30 fs. The molecules are initially at T = 8 K. The total 
interaction strength is Ptot = 5. The results are shown as a 
function of the pulse train period r. (a) Absorbed rotational 
energy for a non-rotating (S = 0) pulse train, (b) Projection 
of the molecular angular momentum on Z for a rotating (S — 
7r/4) pulse train. The angular momentum is given in units 
of h. The shaded areas mark time-delays for which the two 
isotopologues rotate in opposite direction after the excitation. 



higher numerical complexity of the problem, we only con- 
sidered delta-pulses. Comparison with experiments (37j 
shows that this approximation is well justified. 

Unlike molecular nitrogen, oxygen has a non-zero total 
electronic spin in its ground state. There is a coupling 
between the electronic spin and the orbital angular mo- 
mentum, leading to splitting of the rotational levels as 
shown in Figure [3J Also, the orbital angular momentum 
A" is not identical to the total angular momentum J any 
more, but J = N - 1, N, N + 1. Note that due to the 
symmetry of the molecule, only odd values are permitted 
for N. 

For oxygen, we define the population Q(N) of a rota- 
tional level A" as 



N+l 



Q(JV)=J> E E i<vW 



(36) 



J=N-1 M=-J 



Here, i denotes the initial state and <?i is the corre- 
sponding statistical weight. The population of counter- 
clockwise rotating states Ql{N) and clockwise rotating 
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FIG. 11. (Colour online) Population Q(N) (see Eq. (f36]>) of 
the rotational levels N = 3 and N = 5 for 16 O2 at T = 8K, 
after interacting with a train of delta-pulses with intensity 
envelope given by ((6J (A = 2). The total interaction strength 
is Ptot = 7.5. Note the different scales for the ordinates. 



states Qr(N) is given as 



E 

. J,M>0 



i(#) = I> E 1 

c%,jnm\ +^/2\ci, 



JN0\ 



, J,M<0 



(37a) 



(37b) 



The directionality of the states with given A^ is defined 

as 



e(N) 



h(N)-Q R (N) 
h(N)+Q R (N) 



(38) 



In Figure [TT] we show the final population of the rota- 
tional levels A^ = 3 and A^ = 5 after the interaction with 
the pulse train. In Figure [12] the directionality of these 
levels is shown. We can see the same basic line struc- 
ture as before for the nitrogen molecules. However, due 
to the different pulse train and because of the level split- 
ting, the lines are broader, and the pattern becomes more 
complex, especially for larger values of r. On the other 
hand, the chessboard pattern seen in the directionality 
plots for molecular nitrogen is not found for molecular 
oxygen. This is due to the fact that the chessboard pat- 
tern is caused by very weak side bands, which do not exist 
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FIG. 12. (Colour online) Directionality e(N) (see Eq. (|38]ll of 
the rotational levels N = 3 and N = 5 for 16 O2 at T = 8K, 
after interacting with a train of delta-pulses with intensity 
envelope given by ((6J (A = 2). The total interaction strength 
is Ptot = 7.5. Note the different scales for the ordinates. 

for molecular oxygen due to the level splitting. One can 
also see that the plots for the level N = 3 look more com- 
plex than the plots for the level N = 5. This is caused 
by the relatively stronger splitting for lower rotational 
levels. Our calculated results very well resemble the ones 
from experiment [37| . 

IV. CONCLUSIONS 

In this paper, we provided a detailed theoretical anal- 
ysis of molecular rotational excitation by a "chiral pulse 
train" , which was introduced in [37| and presented in 
more detail in the accompanying experimental paper [3|| . 
The chiral pulse train is formed by linearly polarised 
pulses, uniformly separated in time by a time-delay r, 
with a constant pulse-to-pulse angular shift S of the po- 
larisation direction. We showed that for certain combi- 
nations of t and 6, molecular rotation with a strong pref- 



erential rotational sense (clockwise or counter-clockwise) 
can be excited. In two-dimensional plots of the excited 
population and the rotational directionality as a function 
of the pulse train period and the pulse-to-pulse polarisa- 
tion shift, a distinct pattern is found. It is made out of 
diagonal lines along which a strong preferential rotational 
sense is achieved. Our analysis shows that this pattern 
is caused by quantum interferences of different excitation 
pathways, which interfere constructively along the above 
mentioned lines. 

We demonstrated the feasibility for selective excitation 
of nuclear spin isomers and isotopologues in a mixture by 
the chiral pulse train. We demonstrated the selectivity 
using para-nitrogen and ortho-nitrogen as an example. 
By choosing the parameters of the chiral pulse train such 
that they address only the states of certain parity, one 
can selectively excite one of the isomers. Since for the chi- 
ral pulse train one can also influence the direction of the 
molecular rotation, it is even possible to induce counter- 
rotation of different nuclear spin isomers. Selective ex- 
citation of isotopologues can be reached by making use 
of the different rotational time-scales of different isotopo- 
logues. The pulse train parameters can be chosen such 
that they lead to strong excitation of a preferable isotopo- 
logue. For other isotopologues in the mixture, the same 
pulse train most likely leads to a destructive interference 
of different excitation pathways, so these isotopologues 
are at best only weakly excited. Spin isomer and iso- 
topologue selective excitation using the chiral pulse train 
was recently shown in experiments (32l . |38| , demonstrat- 
ing a good agreement with our theoretical analysis. 

Finally, we investigated the excitation of the more com- 
plex 16 02 molecule by the chiral pulse train. For this 
molecule, the rotational levels are split due to spin-spin 
and spin-orbit interactions. We also used a slightly more 
complex pulse train as employed in experiment |37| . In 
spite of these complications, our main conclusions remain 
valid also for the oxygen molecule. 
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